tha = [-89:90];
phi = [1:360];
tha = [1:180];
fday = 1000;
B = get_pomme(fday,tha,phi);

for i = 1:360,
L(:,i) = abs(Bz(:,i)) == min(abs(Bz(:,i)));
end;  %dip equator




load topo;

for i = 1:100,
   

[x,y,z] = sphere(50);

cla reset
axis square off
props.AmbientStrength = 0.1;
props.DiffuseStrength = 1;
props.SpecularColorReflectance = .5; 
props.SpecularExponent = 20;
props.SpecularStrength = 1;
props.FaceColor= 'texture';
props.EdgeColor = 'none';
props.FaceLighting = 'phong';
props.Cdata = topo;
surface(x,y,z,props);
light('position',[-1 0 1]);
light('position',[-1.5 0.5 -0.5], 'color', [.6 .2 .2]);
view(3);
pause(0.1);
topo = circshift(topo,[0,-1]);
end;